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Abstract 



In absence of explicit solutions of the perturbation equation of a static 
symmetrical homogeneous space-time, the best we can do is to construct 
a gwasi-transformation. In this framework, we solve the perturbation 
equation with initial data and a number of results are derived. 

Far from the horizon of a black hole of even space dimension N, a 

mass-lcss field decays as r'(— + t^) ^ in space-time, where i is a 
harmonic number of the sphere. 

A relation of energy and momentum of a particle with mass in a hyper 
black hole is discovered and a solution to the equation of Klein-Gordon 
in the metric of Schwarzschild-Tangherlini with initial data on the hy- 
persphere is proposed. Also, the Green's function of the Klein-Gordon 
equation in Schwarzschild coordinates is calculated. This function is a 
sum on the harmonic modes of the sphere. The first term is a double 
integration on the spectrum of energy and the momentum of the parti- 
cle. Far from the horizon, the double integration is approximated by an 
integration on a line defined by the relation of energy and momentum of 
a free particle. From here, the potential of Yukawa is derived. Finally, 
the linear perturbation equations arc derived and solved exactly. 

The master equation with initial data of a small perturbation in a 
static symmetrical homogeneous space-time, like a (possibly higher di- 
mensional) Schwarzschild black hole, is studied. A main statement of 
the article is that for each harmonic mode of the horizon there are two 
solutions that behave similarly at large. In the basic mode, the asymp- 
tote of a field (an eigentensor of the Lichnerowicz operator, for example) 

decays at infinity according to the universal law (— + t^) ^ . These 
solutions occur in an integral form. 

The analysis we present is of a small perturbation to the full static 
symmetrical solution. The higher order perturbations will appear in a 
sequel. We determine independently perturbations of the space-time in 
dimension 1-1- AT > 4, where the system of equations can be reduced to 
a master equation — a tensor differential equation. The solutions are 
integral transformations which in some cases reduce to explicit functions. 



1 Introduction 

The manuscript deals with the evolution of fields in background space-times 
with certain symmetry, mainly the higher-dimensional Schwarzschild black 
holes (or general spherically symmetric curved) geometries in a linearized 
approach. We solve the Klein-Gordon equation and the equation for metric 
perturbation introducing a phase space. A number of results are derived. 
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Since it is an important problem, some motivations are described below. 
The reader should probably like to read the original works by Reggae and 
Wheeler and others, for instance: General theory of perturbation analysis 
and motivation [1] [2] [3] [4], Quasi-normal modes [5] [6], and Late-time 
Tails in black hole backgrounds [7] [8] [9] [10] [11]. 

Space-time is described by 

R 

(1) Ricg = -0, 

where 0, R and Ricg are the metric, Ric scalar and Ric tensor of empty 
universe. The Einstein equations (1) can be derived from the Einstein- 
Hilbert action 

(2) ^ / 

where g = det 0. 

A consistent theory of quantum gravity, such as string theory, seems to 
require the existence of higher dimensions, which in order to not contradict- 
ing observational evidence must be compact on small scales. String theory 
has made some important progress in explaining the entropy of certain black 
holes by counting microscopic degrees of freedom, and thus scenarios with 
extra, compact dimensions must be taken seriously. In turn, some of the 
most interesting objects to be studied within string theory are those that 
possess event horizons, possibly extended in the compact extra dimensions. 

Higher-dimensional theories of gravitation have recently attracted much 
attention, and static symmetrical homogeneous solutions have played an 
important role in revealing various aspects of gravitation peculiar to higher 
dimensional systems [12] [13] [14] [15]. A noteworthy fact is that, com- 
pared to the severely constrained situation in four space-time dimensions 
[16] [17] [18], higher dimensional black holes can have many very different 
varieties [19] [20] [21] [22] [23]. For this reason, it is important to study 
the dynamics of fields [24] [25] [26] [27] [28] and space-time perturbations 
[22] [23] [29] in such black hole backgrounds and to find stable solutions. 
Amongst the various types of static symmetrical homogeneous solutions, a 
natural higher dimensional generalization of a black hole — also known as 
a Schwarzschild-Tangherlini black hole [30] — has been assumed to be (or 
believed to be) stable, like its four-dimensional counterpart, and for this 
reason it has been conjectured that this black hole describes a final equilib- 
rium state of evolution of an isolated system, or a certain class of unstable 
static solutions in higher dimensions. However, the question of the asymp- 
tote of a perturbation of a black hole itself is not a trivial matter, and, to 
our knowledge, it has not yet been fully analyzed in higher dimensions. 
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Our manifold will be the mathematical manifold M^~^^ , the first layer 
of the would-be physical space-time of general relativity. The space-time is 
assumed to be 1 -|- dimensional. It has been common practice in recent 
years to denote n as the number of extra dimensions, besides the usual four. 
We are consistent with Ref. [12]. In the mathematical literature about 
topological spaces, it is always implicitly assumed that the entities of the 
set can be distinguished and considered separately (provided the HausdorfF 
conditions are satisfied) , otherwise one could not even talk about point map- 
pings or homeomorphisms. It is well known, however, that the points of a 
homogeneous space cannot have any intrinsic individuality^ . 

For concreteness, we consider as background the Cartesian product of a 
two-dimensional orbit space M"^ and a horizon M^~^ with A^— 1 compact 
dimensions, where 1+N is the total space-time dimension and > 4 (in 
order to avoid space-times with one or less compact spatial dimensions where 
the presence of a massive source is inconsistent with asymptotic flatness). 

Tangherlini's idea was to consider pure gravity in 1-|-A^ dimensions. The 
(l-l-A^)-dimensional action is simply 




where now g is the 1-|-A^ dimensional metric, and R is the corresponding Ric 
scalar. We have the original first two dimensions labeled with coordinates 
where // = 0,1. The A^— 1 remaining dimensions are compacted on a 
hyper sphere and are labeled by the coordinate x £ S'^'^. 

Now we can expand the metric as a spherical harmonics series of the form 

oo 

(4) 5(x,x)=EEs^™'^(^)^™^W- 

/=0 mi 

We find that we get an infinite number of fields in the first two dimensions. 
Modes with mi ^ correspond to massive fields with mass |m/|/vol5''^~^. 
The zero mode corresponds to a mass-less field. As we take volS^"^ to be 
smaller and smaller we see that the mass of the first massive field becomes 
very large. This means that if we compact on a small enough hyper sphere 
we can truncate to mass-less modes in the two-dimensional theory. We 
can only see the extra dimensions by exciting massive modes which are at 
energies beyond our reach. 



^As Hermann Weyl [31] puts it: "There is no distinguishing objective property by 
which one could tell apart one point from all others in a homogeneous space: at this level, 
fixation of a point is only by a demonstrative act as indicated by terms like this and there." 
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Let us now focus on the zero mode ^^^\x), the metric of AA^. We could de- 
fine g so that 

0(0) is the two-dimensional field in the Schwarzschild- Tangher- 
lini metric. In order that our results are more transparent we introduce some 
notation. We denote by r the "radial" coordinate, this is the coordinate 
away from the horizon. In the vicinity of the static homogeneous horizon 
it is useful to introduce "spherical" coordinates x ^ well. We will actually 
introduce the components of the metric in the following way: 

(5) 0(°) = -goo{r)dt dt + gii{r)dr ® dr, g = g^") + r\ 
where 500 and 511 are two derivable functions of the variable r, and 

(6) \ = h{z)ijdz'dz^ 

denotes the metric of the horizon^ , i. e., it is the metric of the Einstein 
space {M^~^, Ij) whose Ric curvature is 

(7) Ric„ = 
S being the Ric scalar of the horizon. 

We have truncated to the mass-less fields. Tangherlini type compactifica- 
tions can be more complicated than simply compacting on a hyper sphere. 
The important thing is that the extra dimensions are small so that we do 
not excite massive modes. We can truncate to mass- less modes and read off 
the effective theory in the first two dimensions. 

We need not restrict ourselves to just two dimension either. In fact, 
higher dimensions have become very fashionable in the last twenty years, 
mainly due to the success of string theory as a possible quantum theory 
of gravity. At the quantum level, boson string theory is only consistent^ 
in twenty-six (!) dimensions, although this figure is reduced to ten when 
we introduce super symmetry. Furthermore, there are five distinct string 
theories which can be viewed as different elements of an embracing new 
theory, M-theory [34] [35] [36]. M-theory lives in eleven dimensions and 
has eleven-dimensional super gravity as its low energy limit. 

Traditionally we achieve the reduction down to the first two dimensions 
using Tangherlini techniques. If we start with a (1 -|- A^)-dimensional theory, 
we compact on a small {N — 1) -dimensional manifold. Different manifolds 
generally give different effective theories in the first two dimensions. The 
one thing all of these manifolds have in common is that they are very small, 
and compact. 



■^Accordingly, we also refer to the metric (5) as a maximally symmetric static homo- 
geneous space-time [22] [23] [32]. Ishibashi and Kodama briefly comment on such a 
generalized static homogeneous space-time in section 4 of [33]. 

Actually, boson string theory contains a tachyon, but we will ignore that here. 
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There is, however, an alternative to TangherUni compactification. This 
is the idea that we hve on something caUed a braneworld, where the extra 
dimensions can be infinite. 

In 1+N—p space-time dimensions, an event horizon can be topologically a 
sphere s'^~^~p, but when extended to p extra dimensions, it could naturally 
have topologies either S^~^ in which case it is a higher dimensional black 
hole, or S^~^~P x RP being then a black p-Brana or a black string in the 
p=l case [12] [30] [37]. The R^ topology of the transverse dimensions can 
be compacted giving flat toroidal topological spaces T^. 

One of the first steps towards understanding these extended higher di- 
mensional solutions is to investigate their classical stability against small 
perturbations. If a solution is unstable, then it most certainly will not be 
found in nature (unless the instability is secular) and the solution looses 
most of its power. Of course the next question one must ask is what is the 
final stable result of such instability. Now, in {1+N — p)-space-time dimen- 
sions, the Schwarzschild-Tangherlini geometry is stable against all kinds of 
perturbations, massive or mass-less [24] [1] [38] [33] [39] [40] [41] [42]. On 
the other hand, quite surprisingly, Gregory and Laflamme [43] showed that 
this is not the case for higher dimensional black branes and black strings. 
These objects are unstable. The kind of black p-branes originally studied in 
[43] were solutions of ten dimensional low energy string theory with metric 
of the form 

(8) ds'^ = dsg^hwa + dx'dxi, 

where ds^^-^^^^ stands for the (l+A^— p)-dimensional Schwarzschild- Tangher- 
lini line element, the are the coordinates of the compact dimensions, and i 
runs from 1 to p. The total dimension of the space-time is In string or 

super gravity theories one takes = 10 or = 11, but for generality 
one can leave it as a free natural parameter. 

The specific metric we shall be interested in is given by a time-independent 
symmetric metric that can be written in Schwarzschild-Tangherlini coordi- 
nates [44] as (5). 

After this general explanation of our background space-times, we give 
in the next section our formalism for gravitational perturbations, and the 
basic strategy of our analysis of the asymptote of a perturbation. At this 
stage, our arguments are not limited to the case of black holes, but are 
equally applicable to any case of maximally symmetric static homogeneous 
space-times with a non- vanishing cosmological constant. 
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2 Solution of the perturbation equation of a 
static symmetrical homogeneous space-time 

The study of symmetric perturbations of the Schwarzschild-like solution is 
reduced to solving a single partial differential equation (9). We attempt to 
solve this perturbation equation with initial data in a general static sym- 
metrical homogeneous space-time. 

Deriving the analytic formula for the solutions is important and useful 
for further study. The computation of gravitational radiation was done only 
numerically in a preceding study, e.g. Ref [42]. 

The first decision to be made is to choose our coordinates and ansatz for 
the analysis. We choose to work with one coordinate patch for the whole 
space-time. In our coordinates (where the metric in the r line is in the 
form gii{r)) the coordinate size of the horizon is an invariant and hence the 
coordinate patch does not change. 

Consider a function / of the variable r. 

We denote by the Laplace operator in the metric Ij. 

Consider a base {-Ee} of eigen- vectors of A^. 

Let us consider a field with configuration space M^+^ with N even and 
global coordinates (j){r, t, x) described by a singular Euler-Lagrange equa- 



The symmetry of the background guarantees that time and spatial variables 
will separate. In other words, if we rewrite equation (9) in r and t coordinates 
we find that the coefficients of the differential equation do not depend on 
time t, and so the problem can be reduced to one dimension by applying a 
formal transform with respect to the time variable t. 

However, the perturbation equation (9) is not invariant under simulta- 
neous rescaling of time and space variables r i— > rk, t kt. This lack of 
invariance reflects the presence of mass/length scale in the model. The solu- 
tion is most conveniently analyzed in terms of the phase space coordinates 
k, uj, e. 



tion. 



(9) 



( 



1 

9oo{r) 




= 0. 
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In absence of explicit solutions of the perturbation equation, the best we 
can do is to construct a solution by performing a g?/asi-transformation from 
coordinate variables {t, r, x) to phase space variables (w, k, e). On the other 
hand, such transformation must have the remarkable property that, in the 
special regimes with uj{r, t, x) ~ k, the variables (r, t, x) and {k, uj, e) form a 
basis of off-shell observables for the gravitational field even if the solution 
of the perturbation equation is not known. 

3pt 

Basic set-up.We denote by Bq the function of Bessel of order o, where 



(10) o=ye+(^-l + y 

If {Ee} is a discreet set (continuous set), then denotes a sum (an inte- 
gral) . 

3pt For the transform of (j) we have 

(11) Fik,u;,e)=TcPir,t,x), 
with the inverse transformation being 

^ poo . , /"OO 

(12) r^^^'^EeJ exp(y^tujj J Bo{rk)F{e, k,uj)dkduj. 

That is to say, we propose a solution of the equation (9) of the form (12) so 
the derivatives of the original problem will enter as multiplicative terms. 

A vector Ef, satisfies 

(13) A^Ee = eEe. 
The functions of Bessel satisfy 

(14) 2S; = Bo-i-Bi+o, 

(15) —Bo{x) = Bo-i{x)+Bi+o{x). 

X 

We substitute (12) into (9). 
We use (13), (14) and (15). 
(9) becomes 

exp(^y^tc^) j F{e,k,uj) 
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+ Bo(rA:) "-^^ + r /(r) + 



9n{r) \ r \ gu{r) 



2 \ lnp-|'^oo(r)\2( 1+2)^ 



(16) -^^^ + ^ii^if- dfcd.. = 0. 

9oo{r)J gii{r) J J 

We suppose that for each e and each lo, 

(17) limkF{e,k,uj)B'^{x)\^=rk = lim kF{e,k,uJ)B'^{x)\^=rk■ 
k^o fc— ►oo 

We integrate by parts kF{e, k, u!)B'g{x)\x=rkdk. 
(16) becomes 



gii{r) gooir)J gu{r) 



N_ 



Y 



(18) ^_Jn;^^ ^^^^ ^ Q_ 

511 (?^) y 



Therefore, 



iV 



F(e, A;,u;) ( +r f{r) + - + '^"^^^ 



9ii(?^) goo{r) J gii{r) 



5ii(?') y 

We define 

(20) .w = v((/M + ^^+^)»M + ^^^f^j. 

Then Equation 19 is as 

Assume the support of the function F(e, cj, A;) hes in a narrow neighborhood 
of the curve 

Cj2 =cj(/fc)2. 
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Then Equation 21 is approximately equal to 

d 

oF(e, oj, k) - k—F{e, to, k) = 0. 
ok 

Therefore, there exists a function G such that 

F{e,uj,k) = k°G{e,uj). 

Therefore, (12) is as 

(22) r^--'^EeJ exp (V^iu;)G(e,cu) J k°Bo{rk)dkdw. 
Set 

t = 0. 

Then (22) is as 

/"OO /"OO 

(23) r^-^^Ee G{e,uj) k°Bo{rk)dkdiu. 

g J — oo J 

Then (22) solves Equation 9 with initial condition (23). 
2.1 Extinction of the field at infinity 

Now, the usefulness of this approach with respect to previous ones is that we 
can derive the late-time behavior of fields in this background (see Ref. [11]) 
with our formalism. The reader must be aware of the relevant literature on 
the subject, and of what has been done in recent years. 

Suppose that 

(24) lim f{r) = and lim 9** = 1, * = 0, 1. 

r— »oo r— >oo 

Then 

(25) lim uj{k) = k. 

r — >oo 

We calculate 

(-2r)°(-r2+t2)-|-"r(i+o) 



(26) exp(y^tk^k°Bo{rk)dk 



/o ^ ^ V-^ 

The previous discussion applies to a class of globally hyperbolic, topolog- 
ically trivial, non-compact (asymptotically flat at spatial infinity) space- 
times of the type of Christodoulou-Klainermann ones [45]. In them we have 
[46] [47] [48] [49]: 



1) The imposition of suitable boundary conditions on the fields and the 
gauge transformations of metric gravity reduces the asymptotic symmetries 
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at spatial infinity to the asymptotic Poincare group. The asymptotic im- 
plementation of Poincare group makes the general-relativistic definition of 
angular momentum and the matching of general relativity with particle 
physics. 

2) The boundary conditions of point 1) require that the leaves of the fo- 
liations associated with the admissible 1-|-A^ splittings of space-time must 
tend to Minkowski space-like hyper-planes asymptotically orthogonal to the 
(l+A'^)- momentum in a direction-independent way. This property is con- 
cretely enforced by using a technique introduced by [50] [51] for the se- 
lection of space-times admitting asymptotically flat (l-l-A^)-coordinates at 
spatial infinity. 

These facts entail that there is an effective evolution in the mathematical 
time t which parametrizes the leaves of the foliation associated with any 
H-A'" splitting. 

The idea behind a new quantization scheme is to arrive directly to the 
physical Hilbert space by quantizing only the observables of the system and 
treating the configuration variables as c-numbers (like time in the time- 
dependent Schrodinger equation; the momenta become derivatives with re- 
spect to the configuration variables, like the energy is replaced by the time 
derivative). In gravity this scheme would make sense only if the configura- 
tion variables are coordinate-independent. There will be as many coupled 
Schrodinger equations as configuration variables (gauge invariant) and the 
wave function will depend on as many times (besides the standard one) as 
configuration variables. Every line in this parameter space will correspond 
to a configuration of the classical theory. 

3 Decay of a mass-less field far from the horizon of a black 



Now consider the evolution of a mass-less scalar field (j) in the background 
described by (5). The evolution is governed by the curved space Klein- 
Gordon equation 



The metric appearing in (27) should describe the geometry referring to both 
the black Br ana and the scalar field, but if we consider that the amplitude 
of (j) is so small that its contribution to the energy content can be neglected, 
then the metric (5) should be a good approximation to ga/s in (27). We shall 
thus work in this perturbation approach. It turns out that it is possible to 



hole 



(27) 
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simplify considerably equation (27) if we separate the horizon variables from 
the radial and the time variables, as is done in fom' dimensions [52] [53]. 
For higher dimensions we follow [25]. 

The evolution of a minimally coupled scalar field (j) is described by the 
mass-less Klein-Gordon equation (16) [54] [55] [56] [57] [58], with 

/ = 0, 

where (16) is a partial differential operator that contains information about 
the initial shape of the wave packet at t=0. The explicit form of the operator 
(16) is the simplest in this general setting. 

In Section 2, the linear perturbation equation is derived and solved. 

Setting 

e = l{l + N -2), 

l-N ^ 

we see that far from the horizon a mass-less field decays as r\—r + t ) ^ 
in time, where / is the harmonic mode of the sphere. This indicates that it 
behaves like ~ ^i-N-2i fixed r. Ref. [11] clarify that the field behaves 
like ~ ii-^-2« for even A^, and ~ t^-^^-'^i for odd A^. 

The situation for the Tangherlini hyper sphere is clear. Near the hyper 
sphere the fluctuations in the field behave in the same way as for flat hyper 
spheres. However, the field does vanish at the point at infinity, and close to 
it the field perturbations start to decrease. Since this point lies far from the 
hyper sphere we might yet believe that gravity is localized at low enough 
energies. At finite temperature we could hide the point at infinity behind 
another black hole horizon. 

4 The master equation of a homogeneous space-time 

Why a fluctuation about a black hole is regular at infinity? This is an impor- 
tant problem in gravitation and has been studied for long time. Although 
we do not have a definite answer to this question yet, we have rather well 
established answers to the questions of stable (unstable) Gaussian fluctu- 
ations about our black hole, which are summarized as the instabilities of 
gravitation. 

In local particle physics problems in the present universe, the assump- 
tion of the asymptotic flatness will be a good approximation. However, the 
success of the inflationary universe model implies that the cosmological con- 
stant cannot be neglected even locally in the early universe. Further, recent 
developments in unifying theories suggest the possibility that the universe 
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has dimensions higher than four on microscopic scales [59] [60] and that 
mini black holes might be produced in elementary particle processes in col- 
liders as well as in high energy cosmic shower events. Therefore, it is a quite 
important problem whether the regular character of a fluctuation about a 
black hole also holds in the non- vanishing cosmological constant case and /or 
in higher dimensions. 

Here we explain the derivation of the master equation for the perturba- 
tion of a higher-dimensional Schwarzschild black hole. Our purpose is to 
study and present an analytic procedure to obtain solutions for the mas- 
ter equation of a small perturbation in a static symmetrical homogeneous 
space-time. We investigate the asymptote of a linear perturbation of higher 
dimensional space-times in the framework of a gauge-invariant formalism 
for gravitational perturbations of maximally symmetric static homogeneous 
space-times. This formalism was recently developed by Ishibashi and Ko- 
dama [33]. 

In the four-dimensional case, the perturbation analysis of a black hole was 
first carried out by Reggae and Wheeler [1], and its stability was essentially 
confirmed by Vishveshwara [2], followed by further detailed studies [7] [61] 
[62]. The basic observation of those works is that unstable perturbations, 
whose growth in time is unbounded, are initially divergent at the event 
horizon (provided that they are vanishing at large distances) and therefore 
are physically unacceptable. 

A key step in the analysis is the reduction of the perturbation equations to 
a simple, tractable form. For this purpose, in the case of a four-dimensional 
black hole, the symmetry of the background space-time is utilized to classify 
the black hole perturbations into two types, axial (odd) modes and polar 
(even) modes, according to their behavior under the parity transformation 
on the two-sphere. To linear order, the perturbation equations for these 
two types of perturbations are independent, and therefore they reduce to 
a set of two Schrodinger-type second-order ordinary differential equations, 
the Reggae- Wheeler equation for axial perturbations and the Zerilli equation 
[63] [3] for polar perturbations. 

In the higher dimensional case, the role of the two-sphere is played by an 
(A^— l)-variety: According to their tensorial behavior on the (A^— l)-variety, 
perturbations are classified into three types, those of tensor, vector, and 
scalar modes, according to their tensorial behavior on the section of the 
background space-time. The first of these is a novel type of perturbation 
that exists only in the higher dimensional case, while the vector- and scalar- 
type modes correspond, respectively, to the axial and polar modes in the 
four-dimensional case. 
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In order to investigate the asymptote of a perturbation of higher dimen- 
sional space-times, one must study the asymptote of ah the types of pertur- 
bations mentioned above. However, besides from Ishibashi and Kodama's 
work, all analysis made to this time has been carried out for only tensor per- 
turbations [22], and the situation for vector and scalar perturbations has 
not yet been investigated: The asymptote of a perturbation of space-times 
in higher dimensions is not yet established. The treatment of scalar pertur- 
bations^ is the most difficult part of this analysis, because the linearized 
Einstein equations for a scalar perturbation are much more complicated than 
those for a tensor perturbation and do not reduce to the form of a single 
second-order Schrodinger-type in the higher dimensional case. 

Recently, this reduction was carried out by Ishibashi and Kodama [38] in 
the background of maximally symmetric black holes with a non-vanishing 
cosmological constant in an arbitrary number of space-time dimensions. 
From the result of that work, along with the previously obtained ordinary 
differential equations for tensor and vector perturbations [64], we have the 
full set of equations for the perturbations of higher dimensional maximally 
symmetric black holes in the form of a single self-ad-joint second-order ordi- 
nary differential equation — which we refer to as the master equation — for 
each tensorial type of perturbation. These master equations of course, can 
be applied to the perturbation analysis of higher dimensional space-times. 
Let us mention, furthermore, that these master equations have been derived 
in the framework of the gauge-invariant formalism for the perturbations, 
and for this reason, they do not involve the problem of the choice of gauge. 

In the gauge invariant approach we can isolate the configuration variables, 
which carry the descriptive arbitrariness of the theory, from the observables, 
which are gauge invariant quantities providing a coordination of the reduced 
phase space of general relativity, and are subjected to hyperbolic (and there- 
fore causal in the customary sense) evolution equations. In physics the Hole 
Argument is an aspect of the fact that also Einstein's theory is interpreted 
as a gauge theory. 

Recall that we consider as background the Cartesian product of a two- 
dimensional orbit space and a horizon with A^— 1 compact dimensions. 



''The term "scalar perturbation" is often used in reference to a perturbation of a free 
scalar field in works involving quasi-normal mode analysis of black holes. This is concep- 
tually quite different from the scalar perturbations considered in this paper, which are 
modes of space-time perturbations of a vacuum space-time. Furthermore, the effective 
potentials for these two cases are different. Hence, although the quasi-normal mode anal- 
ysis of a test scalar field [24] [26] [27] [28] may shed some light on the problem of the 
asymptote of a perturbation of static homogeneous space-times in higher dimensions, an 
exact analysis of the problem based on the formulation given in this paper is necessary. 
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where is the total space-time dimension and iV > 4 (in order to avoid 
space-times with one or less compact spatial dimensions where the presence 
of a massive source is inconsistent with asymptotic flatness). 

In this background one expects two linearly independent master solutions. 
When the size of the radial variable is small (the order of the size of the hori- 
zon) one expects that the solutions closely resemble a (l+A^)-dimensional 
static homogeneous space-time, while as one increases the distance one ex- 
pects that at some point the perturbed solution will be flat. 

As mentioned above, we have already obtained the master equation for 
each tensorial type of perturbation in the background space-time. Our main 
focus in this section is therefore to prove that for each harmonic mode of 
the horizon there are two solutions that behave similarly at large. In the 

basic mode, this behavior is (— r + t ) ^ . In particular, an eigentensor 
of the Lichnerowicz operator in a (possibly higher dimensional) black hole 
with Zero eigenvalue decays at infinity in this way. These solutions occur 
in an integral form. In particular, there exist no unstable solutions to the 
master equation with physically acceptable initial data. 

Ishibashi and Kodama showed that there exist no well-behaved perturba- 
tions that are regular everywhere outside the event horizon. This is consis- 
tent with the uniqueness [19] of a certain class of black holes in an asymp- 
totically flat background. Although our main interest is in the standard 
higher dimensional black holes, that is, asymptotically flat,^ spherically 
symmetric static vacuum solutions in higher dimensions, we also investi- 
gate the asymptote of a perturbation of other maximally symmetric static 
homogeneous space-times with a non-vanishing cosmological constant with 
respect to tensor- and vector-type perturbations. 

There is a reason to expect a good analytic control of static homogeneous 
space-time perturbations even if we do not have a complete analytic solution 
since we have a good approximation: the space-time is expected to resemble 
closely a (l+A^)-dimensional static homogeneous black hole [30]. 

The motivations for this research are first to obtain a theoretical descrip- 
tion of a small perturbation which is important on its own right, and second 
to gain understanding of perturbation physics through combination with an- 
alytical work. The symbiosis with analytical work comes close to serve as a 
partial substitute of experiments (which are sorely absent in this field): the 



^It has been pointed out, however, that the concept of asymptotic flatness itself involves 
some subtle problems in higher dimensions, especially in the case of an odd number black 
hole dimensions [65]. 
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analytics are essential for understanding static solutions close to perturba- 
tion where the approximation is expected to hold, and serves to formulate 
the aims and methods of the theory. As it turns out, the black holes ob- 
tained numerically show only a single multi-pole mode perturbation to their 
horizon [66] [67], and that lends some hope that the analytic approximation 
would retain some validity for static homogeneous space-times as well. 

Our objective is to describe the perturbation process (to first order). Writ- 
ing down the equations of motion and separating the horizontal variables 
we find that a Ric flatness condition in the horizontal directions yields a 
relation among the radial functions which is similar to a trace condition and 
allows us to simplify the fields. After substitution one can express the fields 
in terms of first and second derivatives. We are left with a second order 
tensor differential equation which reduces to a single second order ordinary 
differential equation in the radial direction, for each metric function (and for 
each harmonic mode) from which the whole metric may be recovered. This 
is the master equation, which after a change of variables simplifies further. 
It would be nice to have a deeper understanding why these reductions were 
to be expected. 

The master equation belongs to the class of Fuchsias equations, where 
Fuchsias means that the equation has only regular singular points on the 
complex sphere which includes infinity, and means that there are exactly 
three such points. As the hyper- geometric case of 3 regular singularities 
there is a general solution to the equation, and several methods are avail- 
able. Our case is a rather special case of the hyper-geometric equation. 
In this case, it turns out that the solutions can be written in terms of an 
integral transformation, and it would be nice to understand why that had 
to be the case. Interestingly, we observe that in some of the relevant cases 
these integral transformations simplify further to explicit functions, and in 
particular, all relevant solutions are explicit functions at infinity (solutions 
which are of even multi-pole number and are regular at the horizon). 

We started working out this problem before we were aware of the results 
of Kodama and Ishibashi [38] [33] [68]. Even though their papers were pub- 
lished more than three years ago, we continued independently after learning 
about them the formalism of gauge invariant perturbation theory, and the 
various changes of variables which are employed there. We were able to do 
so and actually found a somewhat similar master equation. Yet the final 
reduction of our master equation to a hyper-geometric one was motivated 
by those papers. 

In Subsection 4.1 we briefly overview the present status of the investi- 
gations of this problem and discuss what kind of new information can be 
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obtained by the linear perturbation theory. In Subsection 4.2 we discuss 
properties of the eigentensors of the Lichnerowicz Laplacian for a space- 
time with appUcations to the asymptote of a hnear perturbation of higher 
dimensional black holes. In Subsection 4.3 we determine the linearized per- 
turbations of the space-time. We have actually solved the perturbation 
equation with initial data of a static symmetrical homogeneous space-time 
in Section 2. We apply this to perform the analysis of the asymptote of a 
perturbation for tensor, vector and scalar perturbations separately. 



4.1 Higher-dimensional static homogeneous 
space-times 

In this subsection we write the static Einstein equations in a form that will 
be convenient for expansion around a background space-time. This type of 
expansion is known as gauge invariant perturbation and we follow here the 
usual conventions for this type of expansion. 

We denote in this subsection the (l-l-A^)-dimensional space-time metric 
by gij. For the gauge invariant perturbation it is convenient to write the 
Ric tensor in the following form [69] 

1 d'^g^^ 

(28) RlCg = - -gikgjig'^'' dx^dx"^ ^ ^T^'^hmn - Tij , 

where Fj^mn and F^"*^ are the Christoffer symbols of the first and the second 
kind, respectively, and in addition one defines 

_ If dr'- SF^ dg'-^ \ 

(29) r,, ^ -|^5,,_ + 5.,_-g,,g.,__r-j, 

(30) F^- ^ r"r^^n- 

The first step in this procedure is to look at the linearized equations valid 
for weakly gravitating space-times. The metric is taken to be 

(31) 9ij + bij, 

where we denote the perturbation metric by bij. We have (l+A^)(2+A'")/2 
metric functions bij which are functions of (r, t, x)- The linearized field 
equations become 

(32) A^bij + Ricg.^6^' + Ricg^.^&f - 2Riemanng.^^.;6'=' = dj = SirGi+NTij, 
where Ag is the Laplace operator in the metric g. 

We open this subsection with one comment. At higher orders in the per- 
turbation procedure the form of the equations is dominated by the linearized 
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equations and is given by 

(33) Ag<7(-)'*^=F(s,9s), 

where (m) is the order under study and F(g, (9g) are source terms which are 
quadratic, at least, in lower order metric components and their derivatives. 



4.1.1 Generalized static solution 



It is quite easy to find the higher dimensional counterparts of 4D static 
solutions, such as the Tangherlini solution [30]. In reality, it is also possible 
to find a slightly more general family of solutions for the Einstein-Maxwell 
system with cosmological constant by requiring that the space-time metric 
and the electromagnetic tensor can be expressed in the form (5) and 

(34) = ^Eo{eoodt Adt + e^idt Adr + endr A dr). 

The Einstein equations and Maxwell equations determine the two dimen- 
sional metric and the electric field £^0 as [39] [32] 

(35) gu = — , Eq = ^ , 

where 
(36) 

, , S R 2 p""-' ^ ^'Q' 



{N-2){N-1) N{l + N) ' (^_2)(7v_ i)r2(^-2)- 

This formula includes the solution with cosmological constant, which is men- 
tioned above. 

In the special case in which Af^~^ is the sphere S'^~^, the space-time 
becomes static and spherically symmetric. In particular, for 

(37) (^r_2)^iV-i) =^'^ = Q ^ = 0' 
this solution coincides with the Tangherlini solution, and for 

Q 

(38) = 1 and R = 0, 

it gives a higher dimensional counter-part of the Reissner-Nostrum solution. 
More generally, when AI^~^ is a constant curvature space, the solution is 
invariant under SO(iV), ISO(Af-l) and SO(Af-l, 1) for 

P^' (N-2)(N-i) 
respectively. 
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In general, this solution has a naked singularity or does not have a horizon. 
As shown in Ref . [39] , it gives a regular black hole only when the parameters 
Q, R and S are in special regions. In this case, the Einstein space M^~^ 
describes a spatial section of the static homogeneous space-time horizon and 
at the same time the spatial infinity. 



4.1.2 Uniqueness about static solutions 

In the static case with R=0, the uniqueness about asymptotically flat reg- 
ular static solutions are now established in higher dimensions as well: all 
regular solutions are exhausted by the Tangherlini solution for the vacuum 
system [70], Reissner-Nostrum solution given by (5), (34) and (35) with 
(Af_2f(jV-i) ^^'^ R=0 in the non-degenerate case [19] and the higher di- 
mensional Majumdar-Papapetrou solutions in the degenerate case for the 
Einstein-Maxwell system [71]. It is also proved that the Gibbons-Mada so- 
lution and the Tangherlini solution are the only asymptotically flat regular 
static symmetrical homogeneous solutions for the Einstein-Maxwell-Dilator 
system and for the Einstein-Harmonic-Scalar system, respectively [20] [72]. 
However, for R ^ 0, nothing is known about the uniqueness, although it 
was pointed out in Ref. [73] that the approach developed by Anderson may 
also apply to higher dimensional cases with R < 0. 

In the rotating case, the higher dimensional counterpart of the Kerr solu- 
tion is also known. It is given by the Myers-Perry solution [12], which has 
a horizon and is asymptotically flat. In contrast to the static case, however, 
this is not the unique asymptotically flat rotating regular solution in higher 
dimensions, although the uniqueness holds for super-symmetric static ho- 
mogeneous black holes in the 5-dimensional minimal super-gravity model 
[74]. This is because the horizon topology need not be given by the sphere 
in higher dimensions as discussed in Ref. [75]. In fact, Emparan and Reall 
[21] found an asymptotically flat and rotating regular static symmetrical 
homogeneous solution in five dimensions whose black hole surface is home- 
omorphic to x S^. Thus, at present, we have two different families of 
asymptotically flat rotating regular solutions with different horizon topolo- 
gies. Since a more complicated horizon topology is allowed in dimensions 
greater than 5, it is highly probable that other families of regular solutions 
exist. However, near the static limit, it is likely that a uniqueness theorem 
holds. 
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4.2 Fluctuations about static homogeneous space-times 

Fluctuations about static homogeneous space-times relate to their stability. 
In particular, if the cosmic censorship hypothesis does not hold, a regular 
stationary static homogeneous space-time solution will be unstable against 
a generic perturbation, because the formation of a naked singularity must 
be generic [76]. 

In reality, many of the standard static homogeneous space-time solutions 
are known to be perturbatively unstable [77]. 

In contrast to static homogeneous space-times, the fluctuations about 
rotating black holes are established only for the Kerr solutions [78]. We 
do not discuss the fluctuations about these rotating static solutions in this 
article. 

In the proofs of the fluctuations about these static solutions, a key role 
was played by the fact that the basic equations for the perturbations can 
be reduced to a single second-order ordinary differential equation of the 
Schrodinger type, which is often called a master equation. Now, we ex- 
plain how such a master equation is derived and how the asymptote of a 
perturbation is calculated in terms of it. 

4.2.1 Tensorial decomposition of perturbations 

In space-times of dimension the linearizing of the Einstein equations 

yields the following equations for the space-time perturbation bij=6gij: 

(LichnerowicZ|j6)jj — ViVjb + 2V(jV'^6j)fc 

(40) +{-V^V%i + Ag6 + RicJ'6fc;)9u " Y^^^i = '^'^^^^.j, 
where Lichnerowiczg is the Lichnerowicz operator defined by 

(41) (Lichnerowiczg6)ij = A^bij + Ricg-^ft^' + Ricg^^^f - 2Riemanng.^^.;6''^ 

In order to analyze the behavior of the perturbations on the basis of 
these equations, the following two technical problems have to be resolved. 
Firstly, the perturbation variables bij contain unphysical gauge degrees of 
freedom that should be eliminated. Second, these perturbation equations 
are a coupled system of a large number of equations and are quite hard to 
solve directly. 

In the case of our interest, these problems can be resolved with the help 
of the tensorial decomposition and the gauge-invariant formulation [64] [79] 
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[80]. First, the metric of a static homogeneous space-time solution to the 
vacuum or electron-vacuum system in l+A'^ dimensions can be written in 
terms of the two-dimensional orbit space J\f^ as (5). On this background, 
the basic perturbation variables bij can be grouped into three sets, 

• 0< i, j <1, 

• 0< i <1 and 2< j < N+1, and 

• 2< i, j < N+1, 

according to their tensorial transformation behavior on M^~^. Among 
these, bij for 0< i <1 and 2< j < A'^+l, can be further decomposed into the 

gradient part bi and the divergence-free part as 

(42) bi,=Djbi + bff- D,b'^^=0, 

where Dj is the co- variant derivative with respect to the metric hij on M^~^. 
Similarly, bij for 2< i, j < N+1, can be decomposed into the trace part and 
trace-free part as 

(43) bij = bihij + brij-, b^^j = 0' 

and the trace-free part bxij can be further decomposed into the scalar part 
b^\ the divergence- free vector part and the divergence- free and trace- 
free part b^pjj as 

(44) bTij = {Oibj + ^^/j,^. Af,) bf + 2^(,6g ) + b% ; 

(45) Mf^. = 0, b'Jb% = 0. 

By this decomposition, the linearized Einstein equations are decomposed 
into three decoupled sets of equations, each of which contains only variables 
belonging to one of the three sets of variables, the scalar-type variables 
(^00; ^oii bii, bi, bi, 6^^), vector-type variables {bl^\ b^j), and tensor-type 

(2) 

variables {b^lj)- Further, since the co- variant derivative Dj always appears 
in the combination A^=-DiD^ in the perturbation equations, the expan- 
sion coefficients of these variables in terms of harmonic tensors on M^~^, 
which are denoted as (^oo, ^oi, Cii, fP, Hl, H^), {f^ , h!^^) and ^ 
for the scalar-type, vector-type and tensor-type variables, respectively, are 
mutually coupled among only those corresponding to the same harmonic 
tensor (see Ref. [22] [39] for details of this harmonic expansion). Thus, 
the perturbation equations can be reduced to sets of equations on the two- 
dimensional orbit space Af^ with small number of entries independent of the 
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space-time dimension. In addition, after this harmonic expansion, we can 
easily construct a basis for gauge-invariant variables by taking appropriate 
linear combinations of these two-dimensional variables and their derivatives 
with respect to our coordinates of AT^, and the linearized Einstein equations 
can be written as differential equations for these gauge-invariant variables 
[64]. 



4.3 Eigentensors of the Lichnerov^icz operator 

In this subsection we find the linear perturbations for the (l-l-A^)-dimensional 
static symmetrical homogeneous solution and show that an eigentensor of 
the Lichnerowicz operator is regular at infinity. Reggae and Wheeler [1] 
derived the linear equations that describe small perturbations of the four 
dimensional black hole, and here we generalize their method for the static 
case. These perturbations can be interpreted as deviations of the static ho- 
mogeneous black hole from symmetry due to remote masses. In the case of 
a sphere, we are interested in the influence of the static homogeneous black 
hole (images) on itself. 

As explained in the introduction the zeroth order correction in the horizon 
metric is the (l-l-Ar)-dimensional static symmetrical homogeneous black hole 
[30]. More generally, as our background, we consider the (l-l-A^)-dimensional 
metric g of the form (5) on the Cartesian product of the two-dimensional 
orbit space Af^ and the (A^— l)-dimensional maximally symmetric space 
M^~^, with curvature scalar satisfying (■jy_2)^jv-i) "^^^ constant 

scalar of curvature R of the space-time is related to the (l-l-iV)-dimensional 
cosmological constant A as A=2R, and the parameter p^~^/2 defines the 
static solution mass. For example, in the case (jv-2XAf-i) ~^ R=0, this 
mass is given by 

{N-1)p^~^An-i 



(46) 



167rGi+7V 



where ^7V-i=27r^/Vr[iV/2] is the area of a unit (A^— l)-sphere, and Gi+at 
denotes the (l-l-A^)-dimensional Newton constant. We assumed that p > 
to avoid the appearance of a naked singularity in the background. We are 
interested in the static region — called the black hole wedge — in which 

(47) -r^~N N~2 ^ ^ S 



so that t is a space-like coordinate. Thus, when R> 0, the curvature scalar 

{N-2){N~1) 



S must be positive, while when i? < 0, .^_„w^_,^ can be or ±1. The 
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higher dimensional black hole corresponds to the case (jv-2)(Af-i) P > ^ 
and R=0. 

To study perturbations in this background, it is convenient to introduce 
harmonic tensors in the (TV— l)-dimensional symmetric space (M^~^, tf) and 
expand perturbations in terms of them, as we see below. The type of the ten- 
sor defines the type of perturbation. Combining the expansion coefficients, 
one can construct the gauge-invariant variables in the two-dimensional orbit 
space metric spanned by our coordinates y*=(t, r). The linearized Einstein 
equations then reduce to a set of equations for the gauge-invariant quanti- 
ties. It is found from the Einstein equations that there is a master scalar 
variable for each tensorial type of perturbation. Consequently, the linearized 
Einstein equations reduce to a set of three master equations for the master 
variable. Thorough studies of the gauge-invariant formalism for perturba- 
tions and the derivation of the master equation are given in Refs. [38] and 
[64]. 

To determine the eigenvalues of the Lichnerowicz operator we study the 
solutions of 

(48) A^bij + Ric^i^b'^ + Ricg^.;,6f - 2Riemanng.^^.,5'=' = Xbij, 
where b is transverse, normalizable and 

(49) _^ + ^ = 0, 

900 911 



l+N 

(50) h'^bij = 0. 

A variant of this problem has been treated by Reggae and Wheeler [1]. 
They investigated the Lorentz version of this problem with A=0. Their 
methods were refined by subsequent workers, Vishveshwara [2], Zerilli [63], 
Press and Teukolsky [81] [82], Stewart [83], and Chandrasekhar [84]. A=0 
corresponds to a small perturbation of a black hole that remains a classical 
solution. These authors searched for runaway solutions of the form exp(- 
iu;t) X (function of spatial variables with uj complex). They demonstrated 
that lmaginarya;=0 for all solutions of (48) and concluded that black holes 
were classically stable objects. 

We, on the other hand, are interested in solutions to (48), where gij is 
a static symmetrical homogeneous solution and A is not necessarily Zero. 
Positive (negative) values of A correspond to stable (unstable) Gaussian 
fluctuations about our space-time. 
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By applying the above arguments, we obtain three decoupled sets of equa- 
tions for the scalar-type, vector-type and tensor-type perturbations. The 
first two are often called the even modes (or the polar perturbation) and the 
odd modes (or the axial perturbation), respectively. This equation can be 
separated in t, r and M^~^ coordinates, and is exhibited in equations (62), 
(66) and (69). 

We then follow the approach of Reggae and Wheeler [1] and divide the 
space of eigenfunctions. Since the background is static, by the Fourier trans- 
formation with respect to the time coordinate t, 



(51) hj{t, 



Tijit%jir)Eij, ifO<i<l,0<j<l + N,i<j 



where 



(52) 



2 

LO T, 



(53) 



dt2 



"'-ujfjTij, 0<i<l, 0<j<l + N, i<j, 



(54) Af^Eij = eijEij, < i < 1, < j < 1 + TV, i < j, 
and 

(55) AjA, + Ri%/^ + Ric{, ./f - 2Riemannj.^^.^/3'=^ = ^^f5,„ 
for 2< j < iV+1. 



It turns out that the eigenvalues equation of the Lichnerowicz operator sim- 
plify. The linearized vacuum Einstein equations can be brought to a simpli- 
fied form [85]. Because our background space-time is static, with a Killing 
parameter t in the black hole wedge, substitution of these forms into (48), 
together with the conditions of trace- free and transversality applied to bij, 
leads to sets of decoupled ordinary differential equations for the scalar-type, 
vector-type and tensor-type perturbations in the background space-time. 
The corresponding set of equations can be easily casted into a single second- 
order ordinary differential master equation of the Schrodinger type, as first 
shown by Reggae and Wheeler [1]. It was shown by Zerilli that such reduc- 
tion exists for the scalar-type perturbation of a black hole [63]. Although 
these authors derived these master equations by the gauge- fixing method. 
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the master variable can be written in terms of gauge-invariant variables^ 
[38]. Thus, the hnear fluctuations about a black hole are determined by the 
behavior of the functions ^'s for each mode. 

In the subsequent subsections, we examine tensor, vector and scalar per- 
turbations separately, briefly reviewing the space-time perturbations and 
the master equation. 

Lemma 1. The parts of the Riemann tensor are 

Riemanng.^.^; = 

^511 



900 



N-l ( g'u 1 

1 + N r \2gu r 



+900 



R 



1 



l + N 



7r 



911 V^Sii 



(iV-3) 



S 



N - 1 



r 

(56) H {hiihjk - hikhji) + r^Riemanni, 

Lemma 2. We have 

^^"'Sik^j + ^^%kbi - 2Riemanng^^^.;6' 



ijkl 



(57) 
^/ 

(58) 



m + {N-2){N-l)p 
( 2| _ _ 2)p-2+^r-^) h 



25 



2R 

N {N-2){N-iy^ 



+ 2p 



bji — ^Riemanni, h' 

>iikjl 



kl 



0<i,J < 1 

0<i<land2<j<l + N . 
2 < i,j < 1 + N 



Set 

(59) 



1 ^^N-i [900 d ^ 2R ^ 

r^^^y/ 900911 dr \j gii dr N 



^The problem can be characterized by a single dimensionless variable, x—pr ^. 
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4.3.1 Tensor perturbations 

We begin by considering tensor perturbations, which are given by 

..n^ h -/O' 0<i< l,0<i< l + iV,i<i 

^^"^ ^''~\ mj, 2<^,j<l + N 

where the tensor Pij is defined as solution to the eigenvalue problem on the 
(A^— l)-variety, i.e., Ref. (55) with 

(61) Pi = 0, D,pi = 0. 

The eigenvalues ^ are given by l{l+N—2)—2+2{N—l) for (jv_2Xjv-i) ~ 
(the hyper sphere) and form a continuous set of non-negative values for 
S* < 0. Note that, in the four-dimensional black hole (A^=3), there are no 
corresponding modes, as there are no tensors of this type on a 2-sphere. 

It can immediately be seen that the expansion coefficient / is itself gauge 
independent, and it can be taken as the master variable for a tensor pertur- 
bation. The master equation follows from the vacuum Einstein equations 



(62) 



N -2) goo 



4.3.2 Vector perturbations 

Next, consider vector perturbations, which are given by 

..^s h =1 0<i,j <1 or 2 <i,j <l + N 

^ ' \ ^ijEij, < i < 1 and 2 < J < 1 + iV 

where the vector harmonics are introduced as the solutions of 

(64) (A^ - eij)E,j = 0, 

with eigenvalues Cij that are given by + N — 2) with lij > for 

^^^^ {N-2){N-l) = ^' 

and form a continuous set of non-negative values for 5 < 0. 

A gauge-invariant variable can be constructed. Substituting this into 
the Einstein equations, we can obtain the master variable (^ij. 

The radial part of the Einstein equations can then be reduced to 
(66) (l,+ '4-{N- 2)p-2+^r-^ + 4l] ^ = q. 
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4.3.3 Scalar perturbations 

We now examine scalar-type perturbations in the background space-time. 
These perturbations are given by 



(67) bi^ 



" 1 0, 2<j 
where the scalar harmonics Eij are defined by 
(68) (Af, - eij)E,j = 0, 

with the eigenvalues eij given by + N — 2) for ^jy_2)^jv_i) ^^'^ 

forming a continuous set of non-negative values for (jv-2)^Af-i) — ^• 

We find that the linearized Einstein equations can be reduced to the form: 



(69) L;, + ^ + (AT - 2)(iV - l)p-2+^r-^ + ) = 0. 




4.3.4 The field equations 

We have now 2{l+N) radial fields ^(r), iooir)..., Co(i+Ar)(r), iii{r)..., 
'^i(i+Af)(^) defined in (62), (66), (69) for each mode, 

bijdx^dx-^ = (^ooEoodt^ + ^oiEQi{dtdr + drdt) + ^uEudr^ 

(70) + J] iijEij{dx'dx^ +dx^dx')+il3. 

0<i<l,2<j<l+A' 

Recently, Kodama and Ishibashi have shown that for a static charged ho- 
mogeneous black hole in higher dimensions represented by (5), (34) and (35), 
the perturbation equations can be reduced to decoupled second-order master 
equations of the Schrodinger type, as in four dimensions [38] [39]. These 
master equations in higher dimensions, however, have some new features. 
First, there exists no simple relation between vector and scalar perturbations 
like the scalar-vector correspondence in four dimensions, for d=\+N >4. 
This implies that stabilities for scalar and vector perturbations should be 
studied separately. Second, there exist tensor perturbations for d=l+N >4. 

Now, let us see how these new features affect the fluctuations about static 
homogeneous space-times in higher dimensions. 

When M^~^ is a generic Einstein space, the M^~^-coordinates appear 
in the perturbation equations only through the Lichnerowicz operator 

(71) (LichnerowicZ|j/3)jj = Aij/Jjj + Ric^.^fi^j + Ric,j^.^/3f - 2Riemann|j .^^.^/^''^ 
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Hence, by expanding tensor perturbations in terms of the eigentensors (3 of 
Lichnerowiczj , 

(72) Lichnerowiczf,/3 = fi/S; pj = 0, & pij = 0, 

we obtain a single decoupled equation for ^, which can be easily put into 
the form (62). 

4.3.5 Analysis of the asymptote of a perturbation 

We have solved the master equation with initial data of a static symmetrical 
homogeneous space-time (23). 

Consider the function 
(73) 

{N -2){N -l)p-^+^r-^, 0<i,j<l 
-(iV- 2)p"2+^r-^, < i < 1 and 2 < J < 1 + iV , 

2<.,j<l + iV 

then our universal decay law (-r^ + t^)(^~^)/^ of the field behavior at infinity 
in space-time (23) again agrees with Ref. [11] because of Equations (19) and 
(20). 




5 The Green's function of the Klein-Gordon 
equation in Schwarzschild coordinates 

The present section deals essentially with the field theory in curved space- 
time by studying a class of "exact" (integral transforms) solutions to the 
massive Klein-Gordon field equation in the background of a Schwarzschild- 
Tangherlini black hole (the higher-dimensional Schwarzschild solution) in an 
arbitrary number 1 + of space-time dimensions. 

We again use the relations Equations (19) and (20) in this study. They 
are not exact solutions since it is assumed that the amplitude is so small that 
its contribution to the energy content can be neglected. Some justifications 
for this study are presented. 

We also compute and discuss the behavior of the Green's function of 
the Klein-Gordon equation for a free scalar field in the background of a 
curved higher dimensional spherically symmetric Schwarzschild black hole 
with coordinates coupled to the curvature of the background space-time [86] . 
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The Green's function is a sum on the harmonic modes of the sphere. We 
consider the first term in Subsection 5.3. This term is a double integration 
on the energy spectrum and the momentum of the particle. Far from the 
horizon, we can approximate the double integration by an integration on 
a line defined by the relation of energy and momentum of a free particle. 
From here, we derive the Yukawa potential (91) in our formalism. 

Consider the evolution of a massive scalar field (p in the background de- 
scribed by (5). The evolution is governed by the curved space Klein-Gordon 
equation 



where we denote by g the determinant of the metric. The metric appearing 
in (74) should describe the geometry referring to both the black Brana and 
the scalar field, but if we consider that the amplitude of (p is so small that 
its contribution to the energy content can be neglected, then the metric 
(5) should be a good approximation to ga/3 in (74). We shall thus work 
in this perturbation approach. It turns out that it is possible to simplify 
considerably equation (74) if we separate the horizon variables from the 
radial and time variables, as is done in four dimensions [52] [53]. For higher 
dimensions we follow [25]. 

Consider a number ^. The evolution of a minimally coupled scalar field 
<j) is described by the massive Klein-Gordon equation (9) [54] [55] [56] [57] 
[58], with 



where (9) is a partial differential operator that contains information about 
the initial shape of the wave packet at t=0. This formula also includes the 
cosmological constant. The explicit form of the operator (9) is the simplest 
in this general setting. 

The whole section is based on an integral solution of the Klein-Gordon 
equation. In Section 2, the linear perturbation equation is derived and 
solved. Let's find a relation of energy and momentum of a particle with 
mass in a hyper black hole. We use (20) and calculate 



(74) 




(75) 




(2A;2 + m2)r2-^p 



2+N 



+ o{{p/rf~') 



2\/F~+m2 



(76) 




+ o(fe2)+o((p/r)^-i). 
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The first term is the energy of a free particle. 

The second term is the gravitational potential energy between the masses 
m and ^p~'^~^^ . 

5.1 The Green's function 

To solve the inhomogeneous version of the differential equation (9), we must 
construct a Green's function out of the fundamental system of solutions of 
the homogeneous equation 

POO 

(77) Z,ir,t,x) = r^~^y^E,{x)g{e) / re-*-(^'^'^)i?o(rfc)dfe 

Jo 

(78) Z2{r,t,x) = r^--Y.E,{x)g{e) k"^'^^^^'^^^ B,{rk)dk, 

Jo 

where parameters itw of the differential equation depend on k,r,X as given 
by equations (20) and 

f = m^ + R^-X. 
With these functions, the Green's function is constructed as 

Giruh,xi;r2,t2,X2) = (nra)^"^ J]] Se(xi)^:(X2)x 

e 

/oo i-oo i(tiiAj{k,ri,X)—t2Uj{k,r2,X)) 
/ A;2° Bo{rki)Bo{rk2)dkd\. 
-oo Jo ^ 

where the plus-minus sign of lv is taken depending on the arguments. 

We see that F{e, k, u>) (21) already fixes the Green's function (79). 

The Green's function (79) written above is defined on the Schwarzschild- 
like space-time, but could be analytically continued to the whole r < p 
space. 

Once the Green's function has been determined, it is simple to construct 
the later-time evolution of the wave-packet from the initial data using equa- 
tion (79). As the eigen- value problem (9) is not self-ad-joint, the Green's 
function (79) need not be symmetric in its arguments ri, ti, xi and r2, t2, 
X2- Note that the Green's function is calculated from the A-modes and so 
Zi and Z2 depend on A, but we omitted the third argument in Zi(r, t, x) 
for brevity. 
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5.1.1 Schwarzschild coordinates 

In Schwarzschild coordinates, (79) is simphfied. We denote by Pi the Legend 
re polynomial of degree I. 

We claim that the expression 

^ oo 

G{n,ti,xi;r2,t2,X2) = , , V(l + 2/)P;(Xl •X2)x 

V ^ 1=0 

/oo i-co i{tiw{k,ri,X)—t2w{k,r2,X)) 
/ ^ Bi^i{rki)Bi_^i{rk2)dkdX. 
-oo Jo ^ 2 2 

is the Green's function for the Klein-Gordon field in the Schwarzschild space- 
time and utilize it in the subsequent arguments in essential ways. Actually, 
this expression does not have the global time translation invariance, in spite 
of the fact that the background space-time is static. 



5.2 Incoming Wave-packet 



To get back from the frequency w-interval dependence to the physical time 
evolution, one performs the inverse transformation 

;>oo />oo 

(81) (t>{r,t,x) =r^--yZE,{x) e**" / Bo{rk)F{e,k,uj)dkduj. 

J -oo Jo 



The method of the stationary phase deals with the approximate evaluation 
of Fourier-type integrals 

(82) f{t)= F{uj)ex];)[itS{uj)]dio 

J a 

for large positive parameter t. 

Inverse transform integrals (81) are precisely of the above type, with phase 
S{uj)=u> and the large parameter t being the time coordinate. So the sta- 
tionary phase method tells us that the asymptotic t ^ oo behavior of the 
solution 0(r, t, x) is given by singular points of 

r^-^y^Eeix) / Bo{rk)F{e,k,uj)dk 
Jo 

as a function of uj, i. e. the singular point (20). Therefore the study of 
analytic properties of the Green's function (79) play a key role in under- 
standing the late-time evolution of the wave-packet. The possible sources of 
non-analyticity in the Green's function are listed below. 
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• Branch points in the argument of the square root in uj{k, r, A). 

• Pole at A=0 coming from the defining integral. 

The problem with branches of the argument of the square root in uj{k, r, 
A) is absent in the Green's function G when we approximate the double 
integration by an integration on a line defined by the relation of energy and 
momentum of a free particle, that is to say, by an integration on A ~ m'^+k'^. 

Since the r2, t2, X2 space-time is important for the late-time evolution of 
the wave-packet, it is instructive to take a closer look at the approximation to 
the Green's function (79) there. We proceed to apply this into the Green's 
function (79), noticing that the behavior of Zi is fundamentally different 
depending on the sign of A. I /A < 0, the imaginary exponent ituj{k, r, A) 
in (77) and (78) is large and Zi goes to zero for all r, t, x- If A > 0, the real 
part of itio{k, r, A) is negative and hence Zi is finite for r, t, x- Applying 
this into the Green's function and using the asymptotic behavior of Zi, Z2 
near infinity, we obtain, 

G'(ri,ti,xi;r2,t2,X2) = 

(83) ={nr2f--y^E,{xi)E:{x2) 7^ ^Bo{rki)Bo{rk2)dk. 

g Jq ft + m 

The case A > is precisely where the inverse transformation should be in. 

5.3 From field to particle to force 

In field theory we define 

(84) w{J) = -^j j d^xd^yJ{x)G{x-y)J{y), 

where J{x) is the so-called source function. We can choose any J{x) we 
want and by exploiting this freedom of choice, we can extract a remarkable 
amount of physics. 

We can go on to consider some possibilities for J{x) (which we will refer 
to generically as sources), for example, 

J{x) = Ji{x) + J2{x), where Ja{x) = 6'^'^\x,Xa). 

In other words, J{x) is a sum of sources that are time-dependent infinitely 
sharp spikes located at xi and X2 in space-time. (If the reader likes more 
mathematical rigor than is offered here, he is welcome to replace the delta 
function by lumpy functions peaking at Xa- The reader would simply clutter 
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up the formulas without gaining much.) More picturesquely, we are describ- 
ing two massive lumps sitting at xi and X2 on the space-time, which are 
moving [time dependence in J(x)]. 

What do the quantum fluctuations in the field 4>, that is, fluctuations in 
space-time, do to the lumps? If the reader expected an attraction between 
the two lumps, he is quite right. 

W{J) contains four terms. We neglect the "self-interaction" term JiJi 
since this contribution would be present in W regardless of whether J2 is 
present or not. We want to study the interaction between the two "massive 
lumps" represented by Ji and J2. Similarly we neglect J2J2- 

Plugging into (84) and doing the integration over d'^x and d^y we imme- 
diately obtain 

^ 00 

W{J) = - + 2l)Pi{xi ■ X2)x 

^ ^ ^ (=0 

/oo /-oo pi(tia;(fc,ri,A)— t2'i^(fc,»'2,A)) 
/ Bi^i{rki)Bi_^i{rk2)dkdX. 
odJo a 22 

(The factor 2 comes from the two terms J2J1 and JiJ2-) 

Recall that in path integration formalism Z=Ce^^'^'^^ represents e^^ , where 
E is the energy due to the presence of the two sources acting on each other. 
Setting iW=iE we obtain from (85) 

^ 00 

E = + • X2)x 

47r-i/rir2 
V ^ ^ 1=0 

/ k^+^^ ^ Bi i{rki)Bi i{rk2)dkdX. 

00 Jo X 22 

We will see that this energy is negative. The presence of two delta function 
sources, at xi and X2, has lowered the energy. In other words, the two sources 
attract each other by virtue of their coupling to the field (p. We have derived 
a first physical result in quantum field theory! 

We identify E as the potential energy between two sources. Even without 
doing the integrations we see that as the separation xi, X2 between the two 
sources becomes large, the oscillatory exponential cuts off the integration. 
The characteristic distance is the inverse of the characteristic value of k, 
which is m. Thus, we expect the attraction between the two sources to 
decrease rapidly to zero over the distance 1/m. 
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The range of the attractive force generated by the field 4> is determined 
inversely by the mass m of the particle described by the field. 

We consider the first term 
(87) 

47r^/n?i 7-oo7o A 2 2 

This term is a double integration on the energy spectrum and the momentum 

of the particle. 

Far from horizon p ~ 0, we can approximate the double integration by an 
integration on a line defined by the relation of energy and momentum of a 
free particle, that is to say, by an integration on A ~ 

, ^ Poixi • X2) / {m^ + k^^kBi {rki)Bi {rk2)dk. 

We calculate 

1 1 /"°° 
(89) -— r"2Po(x-X2) / {m^ + k^y^kBi{rk) lim r2^Bi{rk2)dk 
An Jq 2 r-2^0 2 

I 2k smirk) , 
90 =-ir TTo \^dk. 

The integral gives 

(91) E = - J-e-"^^ 

The result is as we expected: The potential drops off exponentially over the 
distance scale 1/m. Obviously, dE/dr > 0: The two massive lumps sitting 
on the space-time can lower the energy by getting closer to each other. 

What we have derived is one of the most celebrated results in twentieth- 
century physics. Yukawa proposed that the attraction between nucleons in 
the atomic nucleus is due to their coupling to a field like the 4> field described 
here. The known range of the nuclear force enabled him to predict not only 
the existence of the particle associated with this field, now called the vr 
meson^ or the peon, but its mass as well. As the reader probably knows, 
the peon was eventually discovered with essentially the properties predicted 
by Yukawa. 

In deriving (91) we direct our attention to the region far from the horizon. 
In this case, oo does not depend on r, i.e. lo is reduced to the ordinary relation 



'''The etymology behind this word is quite interesting (A. Zea, Fearful Symmetry: see 
pp. 169 and 335 to learn, among other things, the French objection and the connection 
between meson and illusion) 
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in the case of a flat space-time. This means that our result (91) is correct 
even if gravity is taken into account. 
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